Trapped Fermi Gases in Rotating Optical Lattices: Realization and Detection of the 

Topological Hofstadter Insulator 



R. O. Umucalilari, Hui Zhai^'^, and M. 6. Oktel^] 

1 Department of Physics, Bilkent University, 06800 Ankara, Turkey 
^ Department of Physics, University of California, Berkeley, California, 94720, USA and 
Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, 94720, USA 

(Dated: February 2, 2008) 

We consider a gas of non-interacting spinless fermions in a rotating optical lattice and calculate 
the density profile of the gas in an external confinement potential. The density profile exhibits 
distinct plateaus, which correspond to gaps in the single particle spectrum known as the Hofstadter 
butterfly. The plateaus result from insulating behavior whenever the Fermi energy lies within a gap. 
We discuss the necessary conditions to realize the Hofstadter insulator in a cold atom setup and 
show how the quantized Hall conductance can be measured from density profiles using the Stfeda 
formula. 



Recently, some fundamental models of many-particle 
quantum systems have been experimentally realized us- 
ing trapped ultra-cold fermions. Some of these exper- 
iments, such as transport in optical lattices has 
demonstrated well known effects with improved preci- 
sion, while others, such as spin imbalanced superfluid- 
ity [2, , have provided access to previously unexplored 
regimes. With constantly improving experimental con- 
trol over ultra-cold systems, it is expected that many 
other fundamental ideas can be tested in laboratory for 
the first time. 

A basic problem in quantum mechanics is the dynamics 
of a charged particle moving in a periodic potential under 
a magnetic field. The single particle spectrum depends 
sensitively on the ratio of the flux through a unit cell of 
the lattice to flux quantum. For a tight binding lattice, 
a single band splits into narrow magnetic bands, forming 
a self-similar energy spectrum known as the Hofstadter 
butterfly Q. The gaps in the Hofstadter spectrum form 
continuous regions for a finite range of flux. For a system 
of non-interacting fermions, it was shown by Thouless et. 
al. that whenever the Fermi energy lies in one of these 
gaps, the Hall conductance of the system is quantized 
[5|. This quantization is topological in nature, and the 
quantized Hall conductance is determined uniquely by 
the magnetic translation symmetry 0. This Hofstadter 
insulating phase is a topological insulator that can be 
characterized by the first Chern number. 

Despite its mathematical elegance, the Hofstadter in- 
sulator can hardly be achieved in solid state systems be- 
cause the magnetic field needs to be thousands of Tesla 
in order to create a magnetic flux which can be compara- 
ble to one flux quantum per unit cell [9]. While in some 
experiments super-lattice structures have been used to 
study the splitting of Landau levels under a periodic po- 
tential, the tight-binding regime (strong magnetic field 
limit) has never been experimentally realized [lo| . The 
main purpose of this Letter is to propose an alternative 
way to achieve and to experimentally study this topologi- 
cal insulator by using ultra-cold Fermi gases in a rotating 



optical lattice. We discuss (i) the conditions to realize the 
Hofstadter insulator in rotating optical lattices, and (ii) 
the manifestation of this insulator in real space density 
profile and the method to detect the Hall conductance in 
a cold atom setup. 

The Hamiltonian for a particle in a rotating lattice is 
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where m is the mass of the particle, 17 is the rotation 
frequency, and uj±_ is the transverse trapping frequency, 
p^ = {px,Py), r — {x,y), and Vq is the maximum depth 
of the optical potential created by a laser light with wave 
number k = 27r/A (for counter-propagating laser beams 
lattice constant a is equal to A/2). Photon recoil energy 
En,, defined as h^k^/{2m) is introduced as the energy 
unit in the following discussion. This Hamiltonian can 
be rearranged as 



H = — z — + Vo[sm^{kx) +sm%ky)] 
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When is close to uj±, the third term represents a 
smooth potential in space. The first term describes the 
motion of a particle under a perpendicular magnetic field 
with strength B = 2mc^l/e. This effective magnetic 
field description has led a number of authors to study 
the quantum-Hall type of physics in rotating quantum 
gases 0- Here, we consider the presence of both rota- 
tion and a lattice potential ^ , employing the Hofstadter 
model, which is appropriate for the description of a non- 
interacting gas of fermions in a lattice. 

We now discuss the conditions for simulating the Hof- 
stadter model in an ultra-cold optical lattice system: 

(1) Hofstadter model is expressly single band; the 
motion of fermions in the periodic potential has to be 
well described by a single-band tight-binding model with 
nearest neighbor hopping. This requirement means that 
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there has to be a finite gap between s- and p-bands, and 
the dispersion of the s-band can be well approximated by 
a cosine function. From the band structure calculations 
for optical lattice potentials, one can easily show that this 



requirement is fulfilled when Vq > 3i?R A rotating 
lattice experiment has recently been carried out ■ Al- 
though in this particular experiment Vq < I-Er, which 
is not deep enough to reach the tight binding regime, 
there is no fundamental reason against increasing Vq a 
few more i?R , as it has been routinely done in static lat- 
tice experiments 
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(2) The "magnetic field" has to be strong enough. The 
dimensionless parameter (f> = a^B/^hc/e) is the magnetic 
fiux quantum per plaquette and is connected to the ro- 
tation frequency as, (j) — 2ma^n/h. For rotation to 
create an effective magnetic flux close to one flux quan- 
tum, the rotation frequency fl must be close to -Er/Zi of 
the lattice. In typical optical lattice experiments the re- 
coil frequency is a few kHz, thus, rotation of the lattice 
at hundreds to thousands of Hz would be enough to reach 
the high magnetic field limit. 

(3) To observe the insulating behavior, the tempera- 
ture has to be lower than the gap of the insulator. The 
gap for the Hofstadter insulator is comparable to the hop- 
ping amplitude t. In a moderately deep lattice, for which 
Vb ~ 3 — 7Eji, for instance, t is of the order of 1 — lOnK. 
This is below currently attainable temperatures. How- 
ever, cooling fermions to this regime is not more difficult 
than achieving degenerate Fermi gases in a lattice, and 
the latter is now the major goal pursued in many key 
laboratories in this field. 

(4) The Hofstadter model is a non-interacting one. 
Strong repulsive or attractive interactions can lead to 
either exciton instability or BCS instability of the in- 
sulating phase (which will be discussed elsewhere), and 
therefore diminish its topological behavior. In ultra-cold 
atom experiments a single species of fermions is naturally 
non-interacting due to the Pauli exclusion principle. 

Even when the above conditions are satisfied, there are 
two factors which complicate the correspondence between 
the results of ultra-cold atom experiments and theories 
developed for 'bulk' systems. The first is the presence 
of an external confining potential in all ultra-cold atom 
experiments, i.e. the third term in Eq.([2|). The second 
factor is that transport measurements are usually very 
hard for cold atom systems. Hereafter, we first calculate 
the density profile of a non-interacting Fermi gas in a ro- 
tating optical lattice, in the presence of a smooth external 
potential We find that the presence of the residual 
trapping potential does not preclude the observation of 
the effects of single particle spectrum, as long as it is 
varying smooth enough on the lattice length scale. Sec- 
ondly, we show that one of the most important transport 
properties, namely the Hall conductance which refiects 
the topological nature, can be inferred from the measure- 
ment of the density profile due to the well known Stf eda 
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FIG. 1: Density of states for the Hofstadter butterfly. Darker 
regions imply greater density. Dashed lines represent the tra- 
jectory of local Fermi energy from the center to the edge of the 
cloud, for different values of (p corresponding to those used in 
Fig. [21 namely (j) ^ 1/3, 1/4, 1/7, and 1/10. Regions marked 
by X and ▲ have Hall conductance axy = ±1; and marked by 



o and ■ have 



±2. 



formula [l5|. 

When the residual trapping potential is slowly varying, 
we can utilize the local-density approximation (LDA) in 
which we define a local chemical potential /i; (r) (or Fermi 
energy) as 



fii{r) = fi-V{r), 



(3) 



where in our case V{r) = m(cj| - rj^)H/2. In what 
follows, we shall denote (w^ — fl^) by w^. 

To find the density profile n{r), we simply count the 
number of states below /i;(r) for the corresponding uni- 
form system as 



n{r)^ / deD{e)e{fi,{r)-e). 



(4) 



So, we shall first find the density of states D{e) for non- 
interacting fermions in a periodic lattice under a mag- 
netic field. By using the Landau gauge A = Bxy, one 
can construct the Hamiltonian for this system as follows 



H 



iAi 



(5) 



where a.i (a|) is the fermionic annihilation (creation) op- 
erator at site i and the sum is over nearest neighbor 
sites. Magnetic field affects the Hamiltonian through Aij , 
which is equal to ±27rn0 , if i and j have the same x coor- 
dinate na, and is otherwise, with the sign being deter- 
mined by the hopping direction. The hopping strength 
t can be obtained as a function of Vo /Ej^ from a simple 
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band-structure calculation. The density of states is dis- 
played in Fig. [TJ When is a rational number p/ q, with 
p and q being relatively prime integers, the energy band 
divides into q bands [J]. In this spectrum the gaps form 
continuous regions in <f> — E plane, although the band 
edges are fractal. When fi lies in a gap, the system is an 
insulator, and as one changes 4> and ^, the topological 
nature and the Hall conductance of the insulator do not 
change as long as one remains within the same gap [^. 
The largest two gaps correspond to insulators with Hall 
conductance Gxy = ±1, and the second largest ones have 
Hall conductance axy — ±2 and so on, as marked in Fig. 

m 

To calculate the integral in Eq. ([4]) efficiently, we note 
that if we take q to be sufficiently large, the band-widths 
of each sub-band become negligibly small. Counting the 
number of states then reduces to counting these bands in 
certain intervals. In our numerical procedure, we deter- 
mine the sub-bands with their edge values Sedge- So the 
number of states per plaquette in the two dimensional 
case can conveniently be expressed as 

= ^ 51 ~ ^<^d9e)- (6) 

In all of our calculations we took q = 401, which is a 
prime number allowing p to be successive integers. 4> val- 
ues for other small denominators of q are approximated 
by properly choosing p. For instance, = 1/10 is approx- 
imated by 40/401, 1/4 by 100/401 and 1/3 by 134/401. 

We now present the density profiles for several (j) val- 
ues. To make a connection with experiments, we refer to 
the work in Ref. [3] in which '^''K atoms are stored in 
an optical lattice with lattice constant a = 413 nm. We 
take Vb = 5£^r,, which gives t = 0.066i?R. The param- 
eters at hand yield E^^/h = 45.98 kHz and t/H = 3.035 
kHz. With the choice u; ^ 355 Hz, the gas extends over 
approximately 60 lattice sites in the radial direction, so 
that the assumption of LDA is satisfied. In Figs. O and 
[3] we fixed the number of fermions at 5000. 

When the local chemical potential fJ,i{r) lies in one of 
the gaps, we have dn{r)/dii{r) = because of vanish- 
ing compressibility. Hence, as one can see by comparing 
Fig. [T] and Fig. [21 corresponding to the energy gaps 
in the single particle spectrum, there appear plateaus in 
the density profile. The discernible number of plateaus 
is related to the size of the energy gaps. For instance, in 
Fig. El^a), the plateau with n = 1 is the band insulator 
with completely filled band, which is topologically trivial 
and has vanishing Hall conductance. Apart from that, for 
(j) — 1/3, the chemical potential trajectory passes through 
two gap regions which gives two plateaus with n = 0.333 
and n — 0.667 respectively. While for cj) = 1/4, there are 
totally four sub-bands, but two of them touch at /i = 0, 
so there are also two gap regions corresponding to two 
plateaus with n = 0.25 and n = 0.75. In Fig. [D^b) we 
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FIG. 2: (a) Density profiles for 5000 fermions with (f> — 1/4, 
!^ = 7.2992 kHz, lu± = 7.3078 kHz (solid line) and = 1/3, 
Q. = 9.7809 kHz, uj± = 9.7873 kHz (dashed line), (b) Density 
profiles for 5000 fermions with (j) = 1/10, Q = 2.9197 kHz, 
uj± = 2.9412 kHz (solid line) and <p = 1/7, Q = 4.1605 kHz, 
uj± — 4.1756 kHz (dashed line). Length is measured in units 
of lattice constant a. Density is given in units of particles per 
lattice site. 

choose two 0's with larger q, where there are more gaps 
in the spectrum and therefore more density plateaus. Ex- 
perimentally, the smaller gap one wants to find, the more 
difficult it is, because it requires larger system size and 
lower temperature. 

In Fig. [3] we show the temperature effect on the visi- 
bility of plateaus. We implement the effect of finite tem- 
perature by incorporating the Fermi-Dirac distribution 
into our calculations as 

n2D(wW,r) = V 7J \ ^1 , 1 -(7) 

2q ^ exp[(eedge - \ii{r))lkBT\^\ 

We observe from Fig. [3] that plateaus will be smeared 
out when fc^T > 0.5i. 

As shown by Thouless et. ai, the topological distinc- 
tion of those insulators manifests itself in the Hall con- 
ductance, which should be quantized in units of e'^/h 
Here, we propose a method to read out the information 
about the Hall conductance quantization from the ob- 
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FIG. 3: (Color online) Density profile for 5000 fermions at 
several temperatures when <j) = 1/4. Plateaus become indis- 
cernible when ksT ~ 0.5i. 



tained density profiles. One can observe how the particle 
density within a particular plateau (labelled by the corre- 
sponding gap) changes with the rotation frequency, and 
the change of plateau density with respect to rotation 
frequency is directly related to the Hall conductance of 
that plateau as follows. 

Stfeda obtained a formula for the Hall conductance of 
a two dimensional charged system as 



'xy 



dN 



(8) 



which is valid when the Fermi energy lies in a gap [l5|. 
Here N is the number of levels below the Fermi energy. 
For a neutral gas we can define a similar response func- 
tion = Jx/Fy, where Jx is the mass current in x 
direction induced by a force Fy in the y direction. Then 
the Sti'eda formula for the rotating system can be written 
as 



1 dN 



(9) 



Expressing the rotation frequency in terms of and N in 
terms of the density per plaquette n{r) (as in our density 
plots), one obtains 



m dn 



(10) 



To measure Hall conductance, we first choose two (p 
values, and identify the plateaus in both density profiles 
that correspond to the same gap. The density difference 
of those two plateaus divided by the difference between 
(j) values gives the Hall conductance. We use Fig. [2] as an 
example to show how this procedure works. The plateaus 
are marked by the same symbol as their corresponding 
gaps in Fig. [T] Two plateaus marked by the same symbol 
are the same insulating phase. In Fig. [H[a), for that 
marked by x, An = 0.333 - 0.25, and An/A0 = 1, for 



that marked by An = 0.667 — 0.75, and An/A<^ = 
— 1. In Fig. m^b), one can get the same quantization 
number for the plateaus marked by x and ▲. In addition, 
there are more plateaus, such as those marked by o and 
■ corresponding to the second largest gaps. For that 
marked by o. An = 0.284 - 0.201 and An/A0 = 2, and 
for that marked by ■, An = 0.716 - 0.802 and Aii/Acp = 
-2. 

In summary, we discussed the experimental conditions 
for simulating Hofstadter model in rotating optical lat- 
tices, such as lattice depth, rotational frequency and tem- 
perature. We calculated the density profile in the pres- 
ence of a smooth residual trapping potential, and showed 
how the density plateaus reflect the insulating behavior 
in a " magnetic field " with incommensurate filling num- 
ber. We also propose a method to measure the Hall con- 
ductance from real space density profiles, without doing 
transport experiments. 
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